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We study the transition between a Coulomb phase and a dimer crystal observed in numerical 
simulations of the three-dimensional classical dimer model, by mapping it to a quantum model 
of bosons in two dimensions. The quantum phase transition that results, from a superfluid to a 
Mott insulator at fractional filling, belongs to a class that cannot be described within the Landau- 
Ginzburg- Wilson paradigm. Using a second mapping, to a dual model of vortices, we show that the 
long-wavelength physics near the transition is described by a U(l) gauge theory with SU(2) matter 
fields. 
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I. INTRODUCTION 

In the presence of strong local constraints, certain clas- 
sical systems exhibit so-called 'Coulomb phases', where 
correlation functions have power-law forms and nontriv- 
ial directional dependence*^ 3 -^ These phases are of 
considerable theoretical interest because of their unusual 
properties, and are also of direct experimental relevance. 
Examples include frustrated magnets such as spin ice^ 
and molecular dimers adsorbed onto surfaces^ 

Coulomb phases stand in contrast to ordered phases 
characterized by a broken symmetry and an associ- 
ated order parameter^ Continuous transitions between 
Coulomb and ordered phases present a novel situation 
where the behavior at the transition cannot be cap- 
tured purely in terms of a Ginzburg-Landau theory of 
the fluctuations of this order parameter. Instead, a 
complete description of the transition requires the long- 
range correlations in the Coulomb phase to be taken into 
account^&i&ii^H 

Classical dimer models^^ are among the simplest pos- 
sible model systems that exhibit Coulomb phases, and 
the discovery of a direct transition into an ordered crys- 
talline phase in the dimer model on the cubic lattice has 
stimulated considerable interest ^ 12 ' 15 i 16 i 17 i 18 The ques- 
tion of whether the transition is continuous or first order 
remains controversial, but it is clear that the correlation 
length at the transition is either divergent or at least 
several orders of magnitude larger than the lattice spac- 
ing. The long-distance properties near the transition can 
therefore be described in terms of a continuum theory in- 
cluding only the relevant degrees of freedom. It has been 
suggeste d 12 i 16 ' 17 that the appropriate description is given 
in terms of a noncompact U(l) gauge theory with SU(2)- 
symmetric matter fields, or noncompact CP 1 (NCCP 1 ). 

In this paper, we analyze the transition in the classi- 
cal dimer model on the cubic lattice by using a mapping 
to an equivalent quantum model in two spatial dimen- 
sions. A brief outline of this mapping and the predic- 



tions that result has been presented previously^ We use 
the standard approach of relating classical statistical me- 
chanics in d dimensions to quantum mechanics in d — 1 
dimensions, which in principle provides an exact identity 
between the partition functions in the two cases. 

By an appropriate choice of the mapping, we represent 
the interacting dimers on the links of a cubic lattice as 
hard-core bosons on the sites of a kagome lattice. The 
Coulomb phase then corresponds to the condensed phase 
of the bosons, and the power-law correlations can be un- 
derstood in terms of the coupling to the phonon mode of 
the superfluid. The thermal transition into the dimer 
crystal is equivalent to a (zero-temperature) quantum 
phase transition from the superfluid to a Mott insulator 
at fractional filling. 

Interestingly, this belongs to a class of unconventional 
quantum phase transitions considered by Balents et al. 19 
In these cases, the phases on the two sides of the transi- 
tion have different order parameters, and a naive appli- 
cation of the Landau- Ginzburg- Wilson (LGW) paradigm 
predicts that a continuous transition requires simultane- 
ous fine-tuning of two independent parameters. Balents 
et alii^ instead propose a critical theory in terms of dual 
vortex degrees of freedom, which allows for a generic con- 
tinuous transition between the two phases. Applying this 
approach to our effective quantum model gives a contin- 
uum gauge theory for the transition that is the same as 
has been obtained using a direct mapping carried out in 
three spatial dimensionsi 16 i 17 

We have previously 11 applied a similar approach to a 
model of nearest-neighbor spin ice, where a transition 
takes place from a high-temperature Coulomb phase to 
a low-temperature saturated phase, which has neither 
power-law correlations nor symmetry breaking^ The 
mapping again leads to a theory of quantum bosons, but 
in that case the thermal phase transition maps to the 
standard quantum phase transition between the super- 
fluid and a vacuum state, described by the conventional 
LGW approach. 

In the remainder of this section, we define the cu- 
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bic dimer model and review its phase structure. In 
Section [TIJ we introduce the mapping from the classi- 
cal dimer model to a quantum model of bosons on the 
kagome lattice. We then show, in Section IIII1 that the 
properties of the Coulomb phase of the classical model 
can be understood in terms of those of the condensed 
phase of the quantum bosons. In Section IIV1 we ad- 
dress the phase transition and use a dual picture in terms 
of vortices to derive a continuum theory to describe the 
critical properties. We conclude in Section[V]with discus- 
sion. In the Appendix, we briefly consider modifications 
to the dimer model that lead to the appearance of an 
intermediate disordered phase and show how this can be 
understood in terms of the quantum mapping. 

A. Model 

We treat a model of classical dimers on the links of 
a simple cubic lattice. In a given configuration of the 
classical model, each link is occupied by either one dimer 
or none, with the close-packing constraint that every site 
of the lattice has a total of precisely one dimer on the 
adjoining links. We define the variables d^r) £ {0,1}, 
giving the number of dimers on the link joining the sites 
at r and r + 8^, where 8^ (/i £ {x, y, z}) is a basis vector 
of the cubic lattice. The close-packing constraint can 
then be expressed as 

^K(r) + ^(r-y=l, (1) 

M 

for all sites r. 

Each configuration is assigned an energy (and hence 
Gibbs-Boltzmann weight) that favors the parallel align- 
ment of the dimers on neighboring links. In the simplest 
case, the energy of a configuration is £ = —Tin, where 
nil is the total number of plaquettes (of any orientation) 
with parallel dimers. In a system with periodic bound- 
ary conditions and an even number of sites in all three 
directions, the minimum of £ occurs when all dimers are 
placed on parallel links, giving n\\ = N , the total number 
of sites. There are six such configurations; one example 
has c£ M (r) = 1 for fi = z and r z odd, and zero otherwise. 
We expect our continuum theory to be equally applica- 
ble to other potentials with the same symmetry that also 
favour columnar crystalline order. 

At temperature T = 0, the system minimizes £ by 
selecting one of these six configurations, breaking both 
the translational and rotational symmetries of the lattice. 
We define the 'magnetization' order parameter 

m^r) = ^(-iy[d^r)-d^r-8 fl )] , (2) 

so that the six ground states have m £ {±8 X , ±8 y , ±8 Z } 
for all r. For small positive T, a low-temperature expan- 
sion predicts that (m) will remain nonzero and directed 
along one of the cubic axes. 



In the opposite limit, T — > oo, there is equal statisti- 
cal weight for all configurations obeying the close-packing 
constraint, the number of which grows exponentially with 
N. This limit has been considered by Huse et al.^ who 
showed that the system exhibits a Coulomb phase, where 
the correlation functions are algebraic at long distances. 
In particular, for the connected part of the dimer corre- 
lation function, one finds the standard 3D dipolar form, 3 

(d,(rK(0)) c ~, r fc^, (3) 

where r/ r = (— l)^f r ^ is ±1 on the two sublattices. This 
form for the correlation functions is expected to persist 
for large finite temperatures. 

The high- and low-temperature phases cannot be 
smoothly connected and must therefore be separated by 
one or more phase transitions. High-precision Monte 
Carlo simulations show that there is in fact a single phase 
transition at a critical temperature^ Tq ~ 1.675, and 
show that this is either continuous or very weakly first 
order. In either case, the correlation length at the tran- 
sition is much larger than the lattice spacing and so a 
continuum description should be expected to capture the 
long-wavelength properties near the transition. 

As we have noted, we expect our theory to be equally 
applicable in the presence of modifications that main- 
tain the symmetry of the configuration energy and the 
ordered states. One can also consider modifications of 
the model that reduce the cubic symmetry of the lattice, 
thereby reducing the symmetry of the effective quantum 
Hamiltonian and changing the degeneracy of the ordered 
states^ We consider one example in the Appendix, in 
which the result is the appearance of an intermediate 
phase between the Coulomb and ordered phases. 



II. MAPPING TO KAGOME BOSONS 

We now describe the first stage of our derivation of 
the continuum theory for the transition, in which we 
map from the statistics of classical dimers to a quantum 
model. We do so by using the standard mapping between 
classical statistical mechanics in d dimensions and quan- 
tum mechanics in d — 1 dimensions, in which one spatial 
dimension of the classical problem is interpreted as the 
(imaginary) time direction for the quantum problem. 



A. Definition of mapping 

In defining the mapping, we have the freedom to choose 
the time direction for the quantum problem, and we do 
so in a way that does not distinguish between the order- 
ing patterns in the low-temperature phase. While one of 
the cubic axes might seem to be a natural choice, tak- 
ing this would necessarily distinguish those cases where 
(m) is parallel to the time direction from those where it 
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FIG. 1: (color online) Projection of the cubic lattice onto a 
(111) plane, with the kagome lattice superimposed. The cu- 
bic sites are shown by the numbers 0, 1 and 2, giving the 
quantum imaginary-time coordinate r mod 3 (see main text). 
Points with solid lines show the sites of the kagome lattice, sit- 
uated at the centres of the cubic bonds (dashed lines) between 
sites with r mod 3 = 1 and 2; they therefore lie in planes with 
r mod 3 — | . The larger red circles shown superimposed on 
some of the kagome sites illustrate the occupied sites in one of 
the six degenerate ordering patterns, corresponding to the six 
ordered states of the cubic dimers. The elementary unit vec- 
tors of the kagome lattice, ei and e2, are shown with dashed 
blue arrows. The coordinate axes shown in the bottom-right 
of the figure are the projection of the cubic axes onto the 
(111) plane. 



is perpendicular. Instead, we choose the [111] direction 
and define the quantum imaginary time as r = r M . 

The mapping follows the standard procedure of using 
a transfer matrix to connect the degrees of freedom in 
one layer of the system to those in the next, followed by 
interpretation of the transfer matrix as the exponential 
of a quantum Hamiltonian. We first divide the links of 
the cubic lattice into layers by the imaginary-time coor- 
dinates of their midpoints. Each layer is treated as a time 
slice and the rows and columns of the transfer matrix T\ 
are labeled by the configurations of two adjacent layers. 
The configurations of a given (111) plane are mapped 
onto the basis states of a quantum Hilbert space by sim- 
ply identifying the presence (or absence) of a dimer with 
the presence (absence) of a boson. 

In order for the resulting quantum model to describe 



lattice bosons, we must ensure two features. First, we 
require the sites of the lattice in one time slice to corre- 
spond to those in the next, and second, we require conser- 
vation of particle number. As far as the first requirement 
is concerned, Figure [1] shows that the midpoints of cu- 
bic links with a given value of r form a kagome lattice, 
and that the lattices formed by adjacent layers are dis- 
placed with respect to each other. Our first requirement 
can nonetheless be satisfied by taking a product of (any 
multiple of) three elementary transfer matrices T\ to de- 
fine a single time step. As illustrated in Figure [TJ planes 
separated by three units in the time direction coincided 
(Each of the three elementary matrices has a different 
form because of the relative displacements between suc- 
cessive kagome layers, but for simplicity we denote them 
all without distinction by T\ . ) 

We also require conservation of particle number, mean- 
ing that the number of bosons in a given time slice should 
be constrained to equal the number in the following. The 
close-packing constraint in Eq. |T]) implies that 



re(r-l) L 



l-£^(r) 



(4) 



where r G r indicates a sum over all cubic sites in 
imaginary-time slice r. This equation states that if the 
kagome plane at r — \ has n bosons and a total of A 
sites, then the plane at r + i will have ^ — n bosons. To 
define a transfer matrix that conserves particle number, 
we must therefore take a product of an even number of 
elementary transfer matrices T\. 

To satisfy the two requirements, we must therefore de- 
fine the full transfer matrix T = T^ T with St an integer 
divisible by both 3 and 2, and so we take St = 6. The 
transfer matrix T therefore has rows and columns la- 
beled by the configurations of two planes separated by 
St = 6, and its elements give the statistical weights for 
those configurations, summed over all possible configura- 
tions of the intermediate planes. 

The partition function for the classical problem is given 
by Z — TrT L / Sr , where L is the length of the system in 
the [111] direction and periodic boundary conditions are 
assumed in this direction. The effective quantum Hamil- 
tonian TL is defined by T — e~ nSr , so that Z is given by 
the quantum partition function at inverse temperature 
[3 oc L. The classical thermodynamic limit is therefore 
given by the quantum zero-temperature limit, (3 — > oo, 
and we will always work in this limit. 



B. Quantum Hamiltonian 

For a finite lattice, it is in principle possible to find the 
transfer matrix T exactly, by considering all allowed con- 
figurations of two planes separated by St = 6 and sum- 
ming the Boltzmann weights of all possible arrangements 
of the intermediate planes. The quantum Hamiltonian ri 
can then be found by taking the (matrix) logarithm of 
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T. For even fairly small lattices, however, the number 
of configurations is large, making this a computationally 
difficult problem, and we have not attempted to find T 
or H exactly. 

Since we are interested in the long-wavelength proper- 
ties of the model near the transition, we will instead use 
general considerations such as symmetry to determine 
the form of the Hamiltonian. As we have noted above, 
the Hamiltonian describes the dynamics of bosons on a 
kagome lattice, with conserved particle number. The 
hard-core nature of the dimers implies that the bosons 
have a similar hard-core constraint, restricting occupa- 
tion numbers to zero or one on any site of the lattice. 

Further, the close-packing constraint in Eq. (p}, which 
ensures particle-number conservation, also implies that 
any triangle of the kagome lattice, of either orientation, 
can be occupied by at most one boson. This implies a 
nearest-neighbour repulsion of infinite magnitude. 

We define the number operator for kagome site i as rij, 
and (hard-core) bosonic creation and annihilation opera- 
tors b\ and b i . The most general quantum Hamiltonian, 
with conserved particle number and obeying the close- 
packing constraint, can then be written in the form 

h = —fj> Uj + niirii -i)+u^2 ninj 

i i (ij) 

where the ellipsis represents other terms, such as three- 
body interactions and correlated hopping terms?!- 1 whose 
precise form is not important. The hard-core constraints 
have been represented by interactions of strength U — > 
oo, the coefficients Vij describe a further- neighbor repul- 
sion, and tij is the hopping. While the Hamiltonian con- 
serves particle number, summing over all classical con- 
figurations means that all particle-number sectors should 
be included in the quantum-mechanical trace. The effec- 
tive quantum problem is therefore defined in the grand- 
canonical ensemble, with a chemical potential fi that, like 
the other coefficients, emerges as an effective parameter. 

A generalization of the dimer models to the case 
where the occupation number of a given link is allowed 
to take on values other than or f would result in a sim- 
ilar quantum Hamiltonian, but with a correspondingly 
expanded on-site Hilbert space. An alternative general- 
ization of the dimer model allows for 'monomers', where 
the close-packing constraint in Eq. {TJ) can be violated to 
permit a site touched by zero or multiple dimers (with 
a finite energy cost). Such a modification breaks con- 
servation of particle number and adds to TL terms such 
as X^W^l + Jfbi) which eliminate the phase-rotation 
symmetry. 

Note that, as is generally the case for effective quantum 
Hamiltonians used to describe classical partition func- 
tions, it is not necessary for TL to be hermitian^ For 
example, the hopping coefficients are in general not sym- 
metric, ^ fa, following from the fact that choosing 



a particular (111) plane to define the quantum problem 
breaks time-reversal symmetry. As we have discussed 
previously in regard to a related problem in spin icejii 
the nonhcrmitian terms are crucial for reproducing the 
correct spatial dependence of the long-range correlation 
functions. 



Locality 

For the analysis that follows, an important condition 
on the Hamiltonian TL is that it should be local, at least 
when projected into a suitable subspace of low-energy 
states. While the classical configuration energy £ and 
the close-packing constraint in Eq. (JTJ are local, this does 
not necessarily imply the same for the transfer matrix or 
quantum Hamiltonian. We have no general proof that 
the locality condition is satisfied, and it is in fact pos- 
sible, by considering states with sufficiently high energy, 
to construct configurations on which the effect of T is to 
cause hopping over arbitrarily long distances. 

We argue, however, that locality is satisfied in the re- 
gion of interest, at low energy near the transition. In this 
region, low-energy configurations of the original cubic 
dimer problem can be described in terms of ordered re- 
gions separated by two-dimensional domain walls, which 
cost an energy proportional to their surface area. The 
intersections of these with a given time slice give the 
one-dimensional domain walls of the quantum problem, 
separating different density- wave orderings of the bosons. 
Consider, in the quantum picture, a time step in which 
a domain wall moves by a large distance, so that one 
of the two neighboring domains grows by an area 8A. 
In the classical partition function, this corresponds to a 
configuration where a domain wall has a section of area 
SA running parallel to the (111) plane, in between the 
two consecutive time slices. Such a configuration has an 
energy cost that grows linearly with 5 A, and hence has 
an exponentially suppressed contribution to the transfer 
matrix. 

Further confirmation of the applicability of a local 
Hamiltonian comes from the analysis of the Coulomb 
phase in Section HTT1 which reproduces the correct power- 
law form of the long-range correlation functions [see 
Eq. (fT4]) ] on the assumption that the low-energy exci- 
tations in the superfluid phase are phonons with a linear 
dispersion. A definitive answer to the question of local- 
ity could of course be found by computing the quantum 
Hamiltonian exactly on a sufficiently large lattice. 

C. Phase structure 

A major advantage of this particular choice of mapping 
is that the six distinct ordering patterns for the classical 
dimers map to six ordered states of the quantum prob- 
lem, related to each other by symmetry. As noted in 
Section II Al the states that minimize the classical config- 
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FIG. 2: (color online) Part of the cubic lattice showing one 
of the ordered states of the dimers (in blue). A (111) plane 
is superimposed on the crystal structure, cutting diagonally 
through the cube. Such planes comprise the 'time slices' on 
which the quantum Hilbert space is defined. The sites of the 
two-dimensional problem are situated where the (111) plane 
intersects links of the cubic lattice; these form a kagome lat- 
tice. Where a cubic link is occupied by a dimer, the corre- 
sponding kagome site is occupied by a boson (shown with red 
spheres) . 

uration energy £ have all dimers parallel to one of the 
three cubic axes, and half of the links of this orientation 
occupied. An example is shown in Figure El along with 
the corresponding arrangement of quantum bosons. 

Cubic links with the same orientation map onto the 
same kagome sublattice, and since half of the cubic links 
of a given orientation are occupied, the same is true of 
the kagome sites of a given sublattice. As illustrated 
in Figure [21 the classical configurations that minimize 
£ map to density-wave states of the bosons, at filling 
(bosons per site) of i. The three states with m = +S Xj y >z 
are related by a rotation of the kagome plane, while m = 
iS z are related by a translation. 

As the (classical) temperature is raised from zero, ther- 
mal fluctuations in the dimer configuration will occur, 
with those of lowest energy being single flipped plaque- 
ttes. In terms of bosons, these correspond to quantum 
fluctuations away from the perfectly ordered density- 
wave patterns and occur once the hopping coefficients 
tij in Eq. ([5]) become nonzero. For sufficiently small hop- 
ping tij relative to the further-neighbor repulsion Vij , the 
ground state remains ordered. We identify this phase, 
where (connected) correlation functions are short-ranged 
and the lattice symmetry is broken, as a Mott insulator 



of bosons with density- wave order. 

When the classical temperature is raised beyond the 
critical value Tq, the dimer order is lost. In the resulting 
high-temperature Coulomb phase, the full lattice symme- 
try is restored and each cubic link has an equal average 
occupation of g. In the quantum model, the transition 
corresponds to a loss of density-wave order at a critical 
hopping tij , and restored lattice symmetry implies a uni- 
form quantum ground state with average particle number 
i on each site. (Note that, since there are three sites in 
the kagome unit cell, this filling corresponds to h per unit 
cell. The hard-core repulsion between nearest-neighbor 
sites also means that filling | is equal to half of the max- 
imum possible filling.) 

This uniform ground state can be identified with the 
superfluid simply by noting that, at fractional filling and 
with neither quenched disorder nor spatial symmetry 
breaking, the only possible phase for quantum bosons at 
zero temperature is a condensate. In Section [ill! we will 
show that the power-law correlations within the Coulomb 
phase are correctly reproduced by the phase mode of the 
condensate, providing further justification for the identi- 
fication of these phases. 

As an aside, we consider the insight into this equiv- 
alence that comes from the phenomenon of off-diagonal 
long-range order (ODLRO). 11 The existence of a nonzero 
superfluid order parameter implies that, in the condensed 
phase, the quantum expectation value ■) approaches 
a nonzero constant in the limit of large separation of the 
points i and j. By contrast, no such ODLRO exists in 
the Mott insulator and the limiting value is zero. 

The quantum expectation value is defined by 

(btb j ) = ^ = ^Tr(T^b\b j ) . (6) 

The quantity Z^j can be understood as the sum over 
histories of the quantum problem, with a particle creation 
event at site i and a particle annihilation at site j, on the 
same (arbitrary) time slice. Returning to the language 
of the classical statistical problem, these events at which 
particle conservation is broken become, according to the 
arguments following Eq. Q , points in three-dimensional 
space where the close-packing constraint is violated. One 
should therefore understand Zij as the partition function 
of the dimer model calculated in the presence of two test 
monomers at positions i and j on the same (arbitrary) 
time slice. (Strictly, one is an empty site, while the other 
is a site where two dimers meet at a site.) 

In the low-temperature phase, such monomers disrupt 
the ordering pattern and cost energy proportional to the 
linear separation between the two sites. By contrast, 
in the Coulomb phase monomers are deconfined, sepa- 
rating them to infinity costs a finite energy^ and Zij 
approaches a nonzero limit for large separation. The 
ODLRO in the quantum superfluid is therefore equiva- 
lent to the deconfinement of monomers; this is consistent 
with our identification of the Coulomb and superfluid 
phases. 
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FIG. 3: (color online) Part of the kagome lattice, shown with 
points joined by solid lines, with the projection of the cubic 
lattice superimposed, as in Figure [1] Four symmetry oper- 
ations, Ki, Kg, R', and Xi, are illustrated in red. The two 
primitive translations Ki and K2 are shown with straight ar- 
rows, while R', a rotation by 60° about the center of a kagome 
hexagon, is shown with a curved arrow. The dashed red ver- 
tical line shows the line of reflection for Xi. The three sites 
of the triangle at the bottom right are labeled a, b, and c to 
denote the three kagome sublattices. 



D. Symmetries of the Hamiltonian 

As we will show in Sections IIIII and HVl an understand- 
ing of the behavior both deep in the Coulomb phase and 
near the ordering transition depends on an analysis of 
the symmetries of the quantum model. We treat here 
the case where the classical configuration energy £ has 
full cubic symmetry. Chen et al.— have studied the ef- 
fects of various modifications that reduce this symmetry, 
and we consider one example in the Appendix. 

First, since the quantum Hamiltonian Tt describes a 
model with conserved particle number, it has a U(l) 
symmetry under global phase rotations of the bosonic 
creation and annihilation operators: bi — > bie li} and 
b\ — ► b\vT l ® . This symmetry is spontaneously broken in 
the supcrfluid phase, where bi has a nonzero expectation 
value. 



Besides this internal symmetry, there are also spatial 
symmetries inherited from those of the classical dimer 
model, but modified in important ways by the particular 
choice of the imaginary-time direction. These symmetries 
can be constructed from combinations of three primitive 
symmetry operations, a translation Ki, a rotation R, and 
a reflection Xi. These operations are illustrated in Fig- 
ure [31 along with the translation IK 2 = MKiR -1 . 

We parametrize the positions of lattice sites in the 
kagome planes by orthogonal coordinates x and y: 

5 = \/^( T v _r *) 

(J) 

y = V2^r z ~ ^r x - hr^j , 

where r M are coordinates referred to the cubic axes, which 
take integer values at the cubic lattice sites. 

The two operators Ki and K2 perform translations by 
the elementary unit vectors ei and e 2 of the kagome lat- 
tice, transforming the coordinates x and y according to 

They can be expressed in terms of pairs of translation 
operators T M for the cubic lattice, chosen so that the 
imaginary time coordinate r is unchanged; for exam- 
ple, Ki = T~ 1 T, / . These transformations map the three 
kagome sublattices to themselves. 

We define R' as a rotation by 60° about the center of a 
kagome hexagon, or, in terms of the cubic lattice, about 
a [111] axis passing though a cubic site with r mod 3 = 
(such as the one at the center of Figure[3|. This maps the 
kagome lattice to itself, but is not a symmetry of the cu- 
bic lattice, since, as can be seen in Figure 02 it exchanges 
cubic sites with r mod 3=1 and r mod 3 = 2. We there- 
fore define the operation R consisting of an improper ro- 
tation by 60° through this axis, which is a symmetry 
of the cubic lattice. In terms of the bosons, R consists 
of the rotation R' followed by a time-reversal operation 
t — ► 3 — t, and it commutes with quantum Hamilto- 
nian TL. (A similar time-reversed symmetry operation 
was found to apply in an effective quantum description 
of spin iceiii The absence of time reversal as an indepen- 
dent symmetry reflects the fact, noted in Section III BI 
that H. is not hermitian.) The rotation R permutes the 
three kagome sublattices cyclically. 

The remaining primitive symmetry operation of the 
kagome model is the reflection Xi through a line running 
perpendicular to the unit vector ei. It transforms x — > 
—x, and exchanges two of the sublattices (labeled a and 
b) while leaving the third unchanged. 

Besides these symmetries of the effective quantum 
Hamiltonian, there are further symmetries of the origi- 
nal classical model that are broken by the explicit choice 
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of the [111] direction as imaginary time. These include 
reflections in the cubic (100), (010), and (001) planes, 
which we denote I x , I y , and I z respectively, and which 
relate different equivalent choices of the imaginary-time 
direction. They cannot be written as operations on the 
quantum Hilbert space, but in terms of the continuum 
space-time action to be derived below, they are simply 
reflections. 

As an aside, one can also consider the operator repre- 
senting translation in the time direction by 3 steps. As 
noted in Section III A| the full transfer matrix T = Tf 
connects two (111) planes separated by 6 steps in the 
imaginary-time direction, and is the simplest choice con- 
sistent with conservation of particle number. We can 
nonetheless consider the operator T 1 / 2 = repre- 
senting 'imaginary-time evolution' by three steps, which 
maps the kagome lattice to itself, and clearly commutes 
with T and hence with the quantum Hamiltonian Ti.. 

Unlike the other symmetries of H, T 1 / 2 cannot be writ- 
ten as a permutation matrix in the occupation-number 
basis, and as for the full transfer matrix, we have not 
attempted to find its precise form. For our purposes, it 
is sufficient to note its effect on the density: as follows 
from the observations of Section [II A[ if a given kagome 
plane has a density of p bosons per site, then the plane 3 
steps later has density | — p. Following our assumption of 
the locality of T, this implies that, in the coarse-grained 
limit, T 1 / 2 simply changes the sign of local density fluc- 
tuations. The microscopic Hamiltonian is not invariant 
under a particle-hole transformation, and this is there- 
fore an emergent symmetry of the long- wavelength limit. 



III. CONTINUUM THEORY FOR COULOMB 
PHASE 

As noted in Section ITAl above the critical temperature 
Tc, the classical dimer model exhibits a Coulomb phase, 
in which there is no ordering, but long-range correlation 
functions have power-law forms and strong spatial de- 
pendence. This behavior can be understood in terms of 
a coarse-grained picture, in which the long-wavelength 
degrees of freedom are described by a solenoidal field. 3 
Deep within the Coulomb phase, this approach predicts 
the dipolar form for the dimer-dimer correlation function 
given in Eq. 

In this section, we will show that the long-distance be- 
havior of the correlation functions can also be obtained 
from the effective quantum model derived in Section HT1 
The power-law behavior of the correlation functions fol- 
lows immediately from the presence of a Goldstone phase 
mode in the superfluid, while the precise spatial depen- 
dence of the dipolar correlations can be reproduced by 
taking into account the symmetries of the effective quan- 
tum Hamiltonian. We have applied a similar analysis to 
a related model of spin ice^i 



A. Kagome continuum action 

To describe the long-distance properties of the super- 
fluid phase of the quantum model, we pass from the mi- 
croscopic description to a continuum action. This action 
is written in terms of bosonic fields '5 corresponding to 
the hard-core boson operator b in the limit where the 
spatial coordinates x and y and the imaginary time r are 
taken as continuous. To preserve the important effects of 
the kagome lattice structure, we define three (c- number) 
fields ^ a corresponding to the three kagome sublattices, 
a € {a, b,c} (illustrated in Figure[3]). 

The continuum action S will contain all powers of the 
fields and their derivatives consistent with the symme- 
tries of the hard-core boson problem, described in Sec- 
tion III Dl We must therefore determine the effects of 
these symmetries on the fields and the coordinates i, y, 
and r. 

Firstly, the symmetry under uniform U(l) phase rota- 
tions of the boson operators leads to the same condition 
on the fields ^f^ and hence restricts the terms in S to 
be those which are invariant under such phase rotations. 
The translation operators Ki and K2 do not affect the 
sublattices and simply lead to translations of the coordi- 
nates x and y. 

Next, consider the operator K, which consists of a ro- 
tation by 60° followed by a time-reversal operation. The 
rotation permutes the three sublattices cyclically and also 
acts on the spatial coordinates, while time reversal cor- 
responds to complex conjugation of the field operators- 1 
^ g — > <!>*. Finally, the reflection Xi exchanges sublat- 
tices a and b while also reflecting the coordinate y — > — y. 

These transformations can be written in a simpler form 



by defining the derivative operators d± 
the vectors d and >&, 



d = 




and 




2 v 



(10) 



and the matrices R and Xi, 



and 




Xi = 




(11) 



The effects of the symmetry operations R and Xi on d, 
and the time derivative d T are summarized in Table [U 

These symmetry considerations allow a continuum ac- 
tion to be written in terms of the field which should 
contain all terms that are invariant under the action of 
the full symmetry group. 

In the Coulomb phase of the classical model, the quan- 
tum bosons condense, so that the field \& acquires a 
nonzero expectation value. This phase breaks no spatial 
symmetries, and so (^a) is equal on the three sublattices 
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TABLE I: Effects of the symmetry operators R and Xi on the 
continuum fields and derivative operators, expressed in terms 
of d and * [Eq. dTOj)] , R and Xi [Eq. {TT}], <j> [Eq. (fT2))], 
and n [Eq. (I13|l ]. as well as the cubic lattice position vec- 
tor r, derivative d, and dimer density field d, discussed in 
Section ITTTBl 
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Rr 


x ir 
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Rd 


Xid 



a. Away from the transition, the long-wavelength prop- 
erties are dominated by the gapless Goldstone mode de- 
scribing uniform rotations of the phase, and correspond- 
ing to the broken U(l) symmetry. Writing ~ e 1 ^, the 
effective action can be expressed in terms of the field 0, 
whose symmetry properties are shown in Table [T] 

The effective phase-only action 5^ can be written in 
the form 



-2 



didydr -<f>(d + <9 2 



(12) 



where the ellipsis denotes terms with higher powers of 4> 
or higher derivatives. (All terms with a single derivative 
can be rewritten as total derivatives and so vanish on 
integration.) This continuum action is explicitly space- 
time symmetric; the relative coefficients of the spatial 
and temporal derivatives are required to be equal [in the 
appropriate units, chosen in Eq. [[7])] by the full cubic 
symmetry of the original model (or equivalently by sym- 
metry under the inversion operators 1^). 

To evaluate the correlation functions of the dimer oc- 
cupation numbers, we must relate these quantities to 
the continuum field ip. First, consider the boson den- 
sity, which we represent by the field n a giving the local 
density measured relative to the average filling of | . The 
symmetry properties of the vector n are included in Ta- 
ble [H note that it is invariant under time reversal. These 
properties are sufficient to identify the density operator 
(to leading order) as 



B. Cubic lattice 

The symmetry operations R and Xi have so far been 
treated as acting within the kagome planes, but they are 
also symmetries of the full cubic lattice. Their action (in 
the continuum limit) on the three-dimensional position 
vector r and derivative operator d are shown in Table [H 
along with the behavior of the field d, which is defined 
as the continuum limit of the dimer occupation number, 
with the average occupation of g subtracted. 

To determine the correlation functions of the dimer 
field d, we must relate it to 4>. It would be consistent 
with the symmetries listed in Table U to identify with 
the combination d^(f>. This is, however, incorrect, as can 
most easily be seen by making use of the cubic reflec- 
tions I M , defined in Section ITlDI changes sign under 
I M , whereas d^ does not. To find the correct relation- 
ship between these fields, consider the effect of the cu- 
bic reflections on the microscopic dimer degrees of free- 
dom, such as d x (0), which gives the occupation num- 
ber for the link between the sites and S x . Under I x , 
this maps to the link between — S x and 0, described by 
dfj,{— <5 x )j in general, the microscopic variable d v (r) maps 
to d v (\ ll Y — 5 lll ,8 ^) under I M . Using i] r , equal to ±1 on the 
two sublattices, we can therefore construct the combina- 
tion 77 r rf At (r), which, after coarse-graining, changes sign 
under I M (since rj r — —rjir—Su for any //). The symme- 
tries in Table U are insufficient in this case because they 
all map from one quantum plane to another; these are 
separated by multiples of St = 6, and so T] T = ( — l) r is 
unchanged. One can instead use the operator T 1 / 2 , de- 
fined at the end of Section III Dl Particle- number conser- 
vation between adjacent planes with t mod 6 = f leads 
to U(l) symmetry under rotations of the phase <j). This 
symmetry therefore acts with the opposite sign in the 
planes with r mod 6 = § + 3, leading to the factor of 

Vr = (-1) T - 

We can therefore identify rj r d^ ~ 9„0, which, together 
with the action given in Eq. (|12[) . allows the Coulomb- 
phase dimer-dimer correlation function to be found. The 
(imaginary-time ordered) propagator for the field cj> is 
simply l/|k| 2 , leading 11 to real-space correlations with a 
dipolar form, 



<d„(r)d v (0)) -Tfc 



|r| 2 <5„ 



(14) 



This form for the correlators was predicted by Huse et 
al.j^ by considering the continuum limit of a coarse- 
grained action for the dimer degrees of freedom. 



d T (, 



(13) 



where the relative coefficient is again fixed using cubic 
symmetry. (The first term is allowed only because of the 
absence of time-reversal symmetry, a consequence of the 
nonhcrmitian nature of ri discussed in Section III Bl ) 



IV. CONTINUUM THEORY OF TRANSITION 

The mapping described in Section |TT] relates the ther- 
mal transition between a dimer crystal and a Coulomb 
phase to the quantum phase transition from a Mott in- 
sulator with density- wave order to a superfluid. A con- 
tinuum theory to describe the phase transition in the 
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dimer model can therefore be found by considering this 
equivalent quantum transition. As noted in Section|H the 
presence of two incompatible order parameters makes the 
standard LGW approach insufficient, and instead, the 
critical theory can be found using a mapping to dual 
vortex fields. 



A. Duality mapping 

This duality mapping has been described in detail by 
Balents et al.,— and we will simply sketch a derivation. 
(Note that Sengupta et ali^I found the critical theory for 
a filling factor of / = -| on the kagome lattice.) The start- 
ing point is a current-loop representation of the quantum 
boson problem,^ where the degrees of freedom are the 
currents J defined on the links of the space-time lattice, 
obeying the continuity equation div J = (where div rep- 
resents the lattice divergence). The essence of the duality 
mapping is a transformation from J to a gauge field A 
on the links of the dual lattice, according to J = curl A. 

It should be noted that mapping to an action written 
in terms of currents, which can be done using a 'Villain 
representation' for the hopping termS ) 19 i 22 involves elimi- 
nating all but the nearest-neighbor hopping. (This can be 
performed explicitly by introducing extra auxiliary fields 
analogous to J to describe further-neighbor processes, be- 
fore integrating these out to give renormalized couplings 
for the currents J.) Reflection symmetries such as Xi 
ensure that the nearest-neighbor hopping coefficients Uj 
are symmetric, and so the nonhermitian nature of TL has 
no effect on the continuum theory near the transition. 
(It was similarly found in a related model for spin ice^> 
that the directed hopping terms in the effective quantum 
Hamiltonian were irrelevant at the transition.) 

In our approach, the bosons occupy the sites of the 
kagome lattice, and so the space-time lattice is not cu- 
bic as in the original dimer problem, but instead consists 
of stacked kagome planes. As illustrated in Figure [4] the 
dual of kagome is the dice lattice ) 21 ' 23 ? 24 and so the gauge 
field Ai is defined on the links £ of a lattice of stacked 
dice planes. The sites of the dice lattice form three sub- 
lattices, which we label as a — 0,1,2, corresponding to 
the three distinct plaquettes of kagome. The positions of 
sites on the dice lattice are given by the (two-component) 
vector x = a\e\ +0262 +cr£, where a\ and 02 are integers 
and C = \{ e \ + e 2) is the displacement between the and 
1 sublattices. 

The currents J, and hence the gauge field A, take in- 
teger values (due to the discrete nature of the quantum 
bosons), and so the model is 'frustrated', in the sense that 
there are many space-time configurations that give nearly 
equal contributions to the action. (This fact is a straight- 
forward consequence of the fractional boson occupation 
number; many arrangements of bosons with filling / = i 
have similar interaction energies.) It is more convenient 
to describe this frustration by introducing 'matter fields' 
ipi on the sites i of the dual lattice, and promote A to 




FIG. 4: (color online) The dice lattice (shown with thin solid 
lines), dual to the kagome lattice (dashed lines). The unit 
vectors ei (horizontal) and e2, as in Figure[2] are shown with 
blue arrows. The thick red lines show two unit cells of the 
dice lattice, or one magnetic unit cell. Within this unit cell, 
the sites of the dice lattice are indicated with black circles 
containing numbers labeling the three sublattices, a = 0, 1, 2. 
The background gauge field At is shown, up to an integer, 
with black arrows, where each arrowhead represents | of a 
flux unit. The arrangement is chosen so that the curl (defined 
as the sum of Ai going counterclockwise around a loop) is 
equal to / = | for every plaquette. (Moving to the right by 
2ei, or one magnetic unit cell, increases At by 1.) 

a continuous-valued field. (This step can be performed 
explicitly by using the Poisson summation formula. 19 ) 
The frustration on A is then lifted, and we shift it by 
a (position-dependent) constant A to make clear that ip 
now carries the frustration. 

The dual theory has a gauge invariance resulting from 
the definition of Ag, and can be written in terms of the 
gauge field Ag and matter fields ipi, 

Sduai = I curlA | 2 
p 

+ $>|^| 2 + ^| 4 ) + --- , (15) 

i 

where J2 P sums over plaquettes p of the dual lattice; 
J2e sums over dual-lattice links £, which start at site i 
and end at site j; and Yli sums over dual-lattice sites. 
Detailed derivations leading to Eq. ([15]) have been given 
by Balents et alj^ 9 - 

The field ipi corresponds to a vortex of the original 
bosons, and the gauge field induces the long-range inter- 
actions between vortices. By duality, the average boson 
density of / = g per site affects the vortices as a flux of 
/ per plaquette. This is represented by the background 
field A obeying curl A = /; we choose the gauge illus- 
trated in Figure |U 

The action 5duai consists of a vortex field ipi with a 
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frustrating hopping term, coupled to a dynamical gauge 
field Ag. Following Balcnts et al.^ our approach will 
be to neglect temporarily the interaction terms, and to 
consider the effect of the frustration on the dispersion of 
a single vortex. The full dual action, including the inter- 
actions, can then be rewritten in terms of the eigenstates 
of the single- vortex Hamiltonian 7i±. The precise form 
of these eigenstates will depend on the details of 6>duai, 
but we will use symmetry considerations to effectively 
block-diagonalize Hi- 



B. Projective symmetry group 

The problem of a single vortex with frustrated hopping 
is equivalent to the Hofstadter problem for a charged par- 
ticle moving on a lattice in the presence of a magnetic 
field. Choosing a specific gauge for the background field 
A reduces the spatial symmetry of the Hamiltonian, and 
it is convenient to introduce a so-called 'projective sym- 
metry group' (PSG) i 19 i 25 We first address the effect of the 
lattice symmetries in the real-space basis, before taking 
the Fourier transform to momentum states, in terms of 
which the eigenstates of Tii can be written. 

The single-particle Hilbert space is spanned by position 
states |x) (where x denotes a lattice site), and Hi is a 
sum of hopping terms of the form 



•We 1 ' 



t(x,b) 



(16) 



for displacements b linking sites of the lattice. 

The PSG associates with a lattice symmetry Q (any 
of the translations, rotations, and reflections defined in 
Section Til D[) . which maps the site x to Q(x), a corre- 
sponding operator Q that commutes with the Hamilto- 
nian Hi. Its action on a state |x) is given by a real-space 
transformation, accompanied by a gauge transformation, 



Q|x; 



= piXQ(x) 



(17) 



Because the Hamiltonian has a lower symmetry than the 
lattice, the operators Q do not obey the group multipli- 
cation table for the full lattice symmetry group. Instead, 
they obey it up to phase factors, determined by the func- 
tions xq- 

Applying Q to the hopping term in Eq. (|16p gives 
Q|x + b)(x|e iQ(x ' b) Q- 1 

= e ±ia(x,b) e i Xe (x+b) e -i X(J (x)|Q ( - x+ b))(Q(x)| , (18) 

where ± is positive (or negative) if Q is an (anti)unitary 
operator. Since Q is a lattice symmetry, there must be a 
corresponding term in Tii given by 

|Q(x + b)) (Q( x )|e iQ (« x )>Q( x + b )-Q( x )) . (19) 

The phases Xt}( x ) should be chosen for all x in order to 
make the expressions in Eqs. (|18p and (fT§]) equal. 



We therefore require 

i[ X Q(x+b)- XQ (x)] = i[a(Q(x),Q(x+b)-Q(x)) T a(x,b)] 



(20) 



This gives a set of equations for xq( x ) which must be 
solved simultaneously. (For the translations K1.2 and ro- 
tation K, a solution can be found for Q a unitary opera- 
tor, but for the reflections such as Xi, Q must be chosen 
antiunitary^) 

Using the choice of gauge illustrated in Figure 0] the 
transformation operators acting on the state |x + cr£), 
with x = a±ei + a 2 e 2 , give 

AT 1 |x + ( rO = (-l) Q2 |x + e 1 + ( 70 (21) 
k 2 \x + aC) = |x + e 2 + CTC) (22) 

R\x + aC) =e T(3«2+5)(2a 1 +« 2 +2i 2 )| t ( x + (j( )j ( 23 ) 

A\|x + aQ) = (-l)4^(«=+ 1 )|X 1 (x + a£)) , (24) 

where R(x+<tC) = (-a2-<5 -i-<5 C r 2 )ei + (ai+a 2 +6 . 2 )e 2 -|- 
aC and X 1 (x + aQ = — (ai + a 2 + <r)ei + a 2 e 2 + crC- Note 
that R exchanges sublattices 1 and 2; we have defined a 
so that = 0, I = 2, and 2 = 1. 

The operators Q form a group with multiplication laws 
equal, up to phase factors, to those of the original space 
group formed by the operators Q. These phases can be 
calculated using the transformations in the position basis, 
giving 



K X K 2 = -K 2 K 1 
K 2 R = c m l 3 RK l 
K 2 R 2 = KiR 2 ^ 



R b 



1 



(25) 
(26) 
(27) 
(28) 
(29) 
(30) 



These commutation properties depend only on the effec- 
tive magnetic flux and are independent of the choice of 
gauge. 

Using the real-space transformations given in 
Eqs. (|2ip - (|2"4"]) . the Fourier transform to momentum 
space can be performed. We define the reciprocal 
lattice vectors e\ and so that • e* = 27r<5y and the 
momentum-space basis by 



aiei + a 2 e 2 + oQ , (31) 



where k = «ie^ + ^e^. The nonuniform phase factors 
Xq (x) cause the operator Q to mix a discrete set of mo- 
mentum values, but at certain high-symmetry points in 
the Brillouin zone (BZ), a smaller set of momenta are 
mixed. 

On the dice lattice with / = i, one finds that a generic 
momentum state belongs to a set of 24 states that are 
mixed, but that there are four points within the (lattice) 
BZ that are closed under the action of the full symmetry 
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group. With our choice of gauge field, their momenta 
are given by k m , n = -(^ + f )e* + + | ) e | with 
m,n S {0,1}, and we expect the global minima of the 
single-particle dispersion to occur at these points^ 

In fact, we argue that two independent linear com- 
binations from the four points form global minima, as 
follows: It can be seen from Figure [4] that K2 remains a 
full symmetry in the presence of the background gauge 
field, and so the corresponding operator K 2 takes a par- 
ticularly simple form. Its action on a position state is 
given in Eq. (|22[) ; acting on a momentum state it gives 
simply 



giving 



K 2 \k,a) =e ie - k |k,a). 



(32) 



This implies that the single-particle Hamiltonian Tli does 
not mix momentum states |k m .„, a) with distinct values 
of n. Considering first n = 0, one of the energy eigen- 
states at the minimum of the single-particle dispersion 
can therefore be written 



|0) = > r„ 



|k m .o, o) 



(33) 



where c ma are coefficients depending on the details of 
Til- (The magnetic BZ is half the size of the lattice BZ, 
so that momenta k m .„ with m = 0, 1 correspond to the 
same point in the magnetic BZ, and are hence mixed by 
-Hi.) 

The state |0) is clearly an eigenstate of K 2 , with eigen- 
value e 17r / 6 . It is straightforward to show, using Eq. ([23]) . 
that the state |1) = e 17r ^ 6 Ki\0) is also an eigenstate of 
K 2 , with eigenvalue — e I7r / 6 , and hence (0|1) = 0. Since 
K\ commutes with TLi (by construction), this state is 
also an energy eigenstate with equal eigenvalue. Using 
the Fourier transform of Eq. (f2"Tj) . one finds explicitly 




(35) 
(36) 
(37) 
(38) 



We have so far considered only the single-particle ki- 
netic terms in the action iSduai given in Eq. (fTS)) . To 
return to the full description, we first define creation op- 
erators v\ and v\ for the single particle states |0) and 
|1). The real-space vortex creation operator projected 
into the low-energy sector can then be written as 



V> f (x) =^(x)4 + <^I(x)u{, 



(39) 



|l)=^c mCT (-l) m |k m , 1 , ( r), 



where ^ n (x) are slowly varying functions (on the lattice 
scale). 

The symmetry properties of the two-component vector 
ip = [ I are determined by those of the states |n), and 

w 

can be summarized as: 



(40) 
(41) 
(42) 
(43) 



Note that the reflection Xi is represented by an antiuni- 
tary operation. 




(34) 



Dual continuum action 



and Ki\l) = e- i7r / 6 |0). 

There are therefore two degenerate minima of the 
single-particle dispersion, \n) with n E {0, 1}. (The same 
result was found by Jiang and Ye^i by directly diago- 
nalizing a Hamiltonian with nearest-neighbor hopping.) 
While the values of the coefficients c ma depend on the 
precise form of the vortex hopping Hamiltonian, the be- 
havior of |0) and 1 1) under the action of the full symmetry 
group can be determined uniquely using Eqs. (125|) - (|30p . 



The action 5d U ai can be rewritten in terms of ip, and 
the other components of the vortex fields integrated out, 
giving renormalized values for the coefficients. The trans- 
formation properties of if reproduce the nontrivial effects 
of the lattice structure and fractional filling /, allowing 
the continuum limit to be taken, with the spatial coordi- 
nates x extended from discrete values describing the dice 
lattice to continuous two-dimensional coordinates. 

The transformations given in Eqs. (|4T?1) (|4"3")) strongly 
constrain the form of the continuum action, which can 
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be expressed as a series in powers of the vortex field ip, 
the gauge field A, and the space- and time-derivative 
operators. As in the Coulomb phase, described in Sec- 
tion Hill the full cubic symmetry of the original dimer 
problem also constrains the final action to be space-time 
symmetric. 

We first consider the terms containing only the vortex 
fields tp. These are most easily found by expressing all 
gauge-invariant (i.e., local and phase-rotation invariant) 
bilinears of the fields <p n in terms of the boson density 
operators, which are the order parameters for the density- 
wave phases. The ordering patterns with which we are 
concerned (such as the one illustrated in Figures Q] and 
[2]) are characterized by nonzero expectation values of the 
momentum-space density at k £ {^e^, ^e^, ^e^ + 

Defining p KlK2 as the boson density operator at mo- 
mentum k = Kiel + the order parameters are P^q, 
p i, and pii. In the continuum limit, we can identify 
these fields with the bilinears of (p n using their symmetry 
properties; we find 



where M = (1 — icr x )/\/2 is a unitary matrix, and <x M 
is a Pauli matrix. In the language of the original dimer 
problem, the order parameter is the magnetization m M , 
related to the dimer occupation number by Eq. ([2]). This 
can similarly be expressed in terms of the vortex fields as 

m M - ¥» t M t er"Mv? . (45) 

We will now identify the allowed interaction terms for 
the vortex fields, and it is convenient to do this by writ- 
ing them in terms of m M . While the full action has cu- 
bic symmetry, the lowest-order terms in fact have SO(3) 
rotation symmetry, which corresponds, via Eq. (|45|) . to 
SU(2) symmetry for the vortex fields (p. Our primary 
concern will be to identify the first term that explicitly 
breaks this higher symmetry. 

Defining the SU(2) Casimir invariant fl = cp^tp, all 
gauge- invariant combinations of the fields <p n can be writ- 
ten in terms of m XiV _ z and f2, and it is straightforward to 
show that |m| 2 ~ J7 2 is also an SU(2) invariant. While 
the action can contain any term involving only f2, terms 
involving functions of ip that break SU(2) are strongly 
constrained by symmetry. 

First note that the only allowed quadratic and quar- 
tic (in the vortex fields) combinations are f2 and il 2 . 
At sixth order, besides fl 3 , one finds the combination 
m x m y m z , which is invariant under all of the operations 
in Eqs. ([4"0)) - P5)) . It is nonetheless excluded by requir- 
ing symmetry under the cubic reflections T M , which take 
"V - * — 7?fy. (In terms of the bosons, this corresponds 
to particle-hole symmetry, which swaps vortices and an- 
tivortices.) 

The lowest-order combination that satisfies all the 
symmetries of the problem, but explicitly breaks the 



SU(2) symmetry, is of eighth order in ip, and its con- 
tribution to the action is 

C 1= vJ2K- ( 46 ) 

Symmetry does not fix the coefficient v, and it is in gen- 
eral allowed to take either sign. In the ordered phases 
that we describe here, however, a single component of 
the magnetization acquires a nonzero expectation value; 
such phases require v < 0. 

D. Emergent SU(2) symmetry 

The full action is given by the continuum limit of 
Eq. (fT5")) and can be written as a three-dimensional inte- 
gral over a Lagrangian density C = Cq + C\ H . The 

lowest order terms Co take the forms expected for a U(l) 
gauge field minimally coupled to matter fields <p, with the 
lattice curl becoming a differential curl, and the modified 
hopping term becoming a covariant derivative: 

C Q = |(V-iA) V | 2 + s | V | 2 + w (|<p| 2 ) 2 + K |VxA| 2 . (47) 

The interaction term C\ is given in Eq. (|46p and the el- 
lipsis represents further terms that are expected to be ir- 
relevant. The cubic symmetry of the original dimer prob- 
lem means that C (considered as the action for a quan- 
tum problem) is space-time symmetric, and this fact has 
been used to express Co in terms of the three-dimensional 
derivative operator V. 

The transition occurs when the quadratic coefficient s 
is tuned through its critical value, and the gauge field 
A acquires a gap by the Anderson-Higgs mechanism. 
The superfluid phase of the boson problem (the Coulomb 
phase of the dimer problem) is equivalent to the Coulomb 
phase of this gauge theory, within which the long-range 
correlations are reproduced by the gapless gauge field. 
The Mott insulator (ordered dimers) corresponds to the 
Higgs phase, where there are no gapless excitations and 
SU(2) symmetry is broken by condensation of the matter 
field ip. 

Note that the gauge field A is by construction noncom- 
pact, so that monopoles are forbidden. Such monopoles 
would give points in space-time with nonzero divergence 
div J of the boson current, which map to monomers in 
the original dimer problem. 

It is not firmly established whether the action C has 
any nontrivial fixed points under the renormalization 
group (RG). As noted by Balents et al.^ such a con- 
jecture is difficult to test analytically, since, for example, 
an expansion in e = 4 — d, where d is the spatial dimen- 
sion, has no weak-coupling fixed points. As we remark in 
Section [Vl the numerical evidence remains inconclusive. 
It is clear from simulations, however, that the transition 
in the original dimer model is, at most, weakly first order, 
with a correlation length of at least hundreds of lattice 
spacings, and so a continuum description is appropriate. 
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The terms appearing in Co have full SU(2) symmetry, 
while the lowest-order term breaking this to the micro- 
scopic cubic symmetry, C\, is of eighth order in the field 
ip. It is therefore highly likely that this term, as well 
as higher-order symmetry-breaking terms, is irrelevant 
in the continuum, and that the effective theory is given 
by Cq and has an emergent SU(2) symmetry. 

This leads naturally to the conclusion that physical 
properties measured sufficiently close to the transition 
should show full SO (3) symmetry, rather than the re- 
duced cubic symmetry of microscopic model. This claim 
is in fact in agreement with qualitative observations made 
by Misguich et al.^ based on their numerical results near 
the transition. First, the measured distribution of the 
order parameter m at the critically becomes increasingly 
spherically symmetric for larger system sizes. Second, 
the dimer correlation function (d M (r)G?„(0)), while tak- 
ing the form given in Eq. ([3]) deep within the Coulomb 
phase, is dominated by a "spinlike" contribution near the 
transition. This follows from the fact that the magneti- 
zation couples directly (without derivatives) to a bilinear 
in the critical field, according to Eq. (|45|) . and forms a 
three-dimensional representation of SU(2). We therefore 
expect 

K(r)m„(0)) ~ S^\r\- d + 2 -^ , (48) 

where the critical exponent rj m is the anomalous dimen- 
sion of the magnetization (and a similar expression holds 
for the dimer correlation function). The absence of a 
weak-coupling fixed point prevents us from making quan- 
titative predictions about the anomalous dimension r\ m . 

As noted by Misguich et al.,— these properties, while 
explained straightforwardly by an SU(2)-symmetric con- 
tinuum theory, are incompatible with other, more obvi- 
ous, candidate continuum theories, such as the 03 model. 

V. DISCUSSION 

In this paper we have presented a derivation of a con- 
tinuum theory to describe the phase transition from a 
dimer crystal to a Coulomb phase observed in simula- 
tions of a classical dimer model on the cubic latticed 
Our approach proceeds by first mapping to a quantum 
model of bosons, which has a corresponding phase tran- 
sition from a Mott insulator with density-wave order to 
a superfluid. The second stage of the derivation consists 
of applying a duality mapping to this quantum model, 
allowing the phase transition to be described in terms of 
an Anderson-Higgs transition for the dual vortex fields. 

The continuum theory we derive coincides with one 
obtained for the same model by Charrier et al.— and by 
Chen et al.,— using a duality mapping applied directly 
to the classical model. The Lagrangian Co is given in 
Eq. ([47]) and is referred to as NCCP 1 ; we expect the long- 
wavelength properties near the transition to be described 
by this theory. In other words, there should be some 
range of length scales, much larger than the lattice scale, 



at which Co provides an appropriate description. This 
assumes that C\ and all higher-order terms are irrelevant 
in the RG sense. As we have noted above, the numerical 
results of Misguich et al. 15 support the claim that at 
least C\ and all other terms breaking SU(2) symmetry 
are irrelevant. 

The question of whether the transition in the dimer 
model is continuous or first order remains open, but is 
clearly related to same question for the transition be- 
tween the Coulomb and Anderson-Higgs phases in the 
model Co- This remains contentious , 26 i 27 i 28 i 29 i 30 i 31 w ith 
numerical evidence so far inconclusive and little prospect 
of insight from analytics. The theory Co has one free 
parameter k (after s is tuned to its critical value), and 
it is possible that it has a continuous transition for some 
value of k and not for others.? (It is not clear how k 
depends on the details of the classical energy £ — for in- 
stance, whether it increases or decreases when one adds 
further-neighbor interactions.) 

In fact, the most economical interpretation of the avail- 
able numerical results is that the original cubic dimer 
model (with only nearest-neighbor interactions) is near a 
tricritical point. Evidence for this comes from the critical 
exponents in the original dimer model and the results of 
adding deformations. 

The critical exponents reported by Alet et al.— for the 
cubic dimer model are, as noted by the authors, con- 
sistent with those expected generically at a tricritical 
point (and in fact those seen at the tricritical point in 
an NCCP 1 model by direct simulations 3 ^). They are in- 
consistent with those for the generic NCCP 1 transition 
as observed in simulations of various models expected to 
be described by this theory . 26 ! 27 ' 28 ' 30 

The scenario that the cubic dimer model with only 
nearest-neighbor interactions is near a tricritical point 
suggests that modifications to the Hamiltonian should 
drive the system away from this tricritical point, cither 
to a strongly first-order transition or to the generic con- 
tinuous transition. Recent numerical results show that 
some perturbations, including ones that preserve the full 
symmetry^ and ones that break it~ can make the tran- 
sition clearly first order. Confirmation of the tricritical- 
point scenario would require the demonstration of a line 
of continuous transitions for some range of parameters. 

Charrier et alt^ have studied a model where Eq. (fl]) is 
still enforced, but the variables (r) are allowed to take 
all integer values (although their simulations treated a 
dual model). As for the unperturbed dimer model, they 
found results consistent with a continuous transition, but 
in this case with exponents in agreement with those re- 
ported for NCCP 1 ! 26 ' 27 ! 28 ' 30 By introducing an energy 
cost for larger values of d M , it is possible to interpolate 
continuously between this generalized version and the 
original dimer model; simulations for a range of values 
would be a useful test of the tricritical-point scenario. 
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FIG. 5: (color online) Left: Illustration of the 'xy' model, in 
which only dimers on the thick blue bonds contribute to the 
energy. The cubic symmetry of the original dimer model is 
reduced and there are only two configurations that minimize 
the energy. Right: Kagome lattice showing the modifications 
to the corresponding quantum model; the sites highlighted in 
blue have an attractive potential energy. The unit cell, shown 
with dashed red lines, is 4 times larger than in the case with 
full symmetry. 
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APPENDIX: MODELS WITH AN 
INTERMEDIATE PHASE 

In their numerical studies, Chen et alJ^ considered var- 
ious modifications of the classical configuration energy £ 
that break the cubic symmetry of the original model and 
reduce the number of degenerate ordered states. While 
certain modifications lead to strongly first-order transi- 
tions, others instead cause the appearance of an interme- 
diate phase between the high-temperature Coulomb and 
low-temperature ordered phases. We will consider one 
particular case, referred to by Chen et as the 'xy' 
model, and show how the phase diagram can be inter- 
preted in terms of the mapping to an effective quantum 
Hamiltonian. 

The 'xy' model is defined so that the only configura- 
tions that minimize the classical energy £ xy are those 
with magnetization m G {d x ,d y }. This can be achieved 
by including interactions only between the subset of 
bonds shown on the left-hand side of Figure[5] Explicitly, 
we can write £ xy = — {n c x + n y ), where 

<= E E^(rK(r + <M ( A -l) 



counts the number of nearest-neighbor parallel dimers on 
even bonds along the [i direction. 

Using the mapping described in Section [Hi the bonds 
included in £ xy correspond to a subset of the sites of the 
kagome lattice shown on the right-hand side of Figure [5] 
In the quantum Hamiltonian 7i xy corresponding to the 
'xy' dimer model, these sites acquire an attractive poten- 
tial energy, reducing the degeneracy of the density-wave 
ordered Mott insulator states that occur in the limit of 
weak hopping. 

Crucially, the reduced symmetry increases the size of 
the unit cell fourfold to 12 sites, so that there are now 
2 bosons per unit cell. It is therefore possible to form a 
zero-temperature insulating state that breaks no symme- 
tries. A simple caricature of such a state has one boson 
localized on every neighboring pair of attractive sites, 
forming a valence bond linking the two. The true ground 
state will involve both fluctuations onto other sites and 
correlations between the positions of nearby bosons. 

Such a valence-bond Mott insulator is expected to be 
the ground state at intermediate ratios of the hopping 
strength to the interactions. When the hopping is in- 
creased, the fluctuations grow larger, leading to a tran- 
sition into a superfluid phase with broken U(l) symme- 
try and gapless Goldstone modes. This transition is de- 
scribed by the standard fixed-density (relativistic) U(l) 
critical theory, as for the Bose-Hubbard model at integer 
filling ) 32 ! 33 in the XY universality class. For smaller hop- 
ping strengths, the further-neighbor interactions between 
bosons will become sufficient to cause long-range order in 
the positions of the bosons. This breaks the symmetry Xi 
(see Figure [3]) giving a density- wave ordered state corre- 
sponding to one of the dimer crystals with magnetization 
m € {S x , d y }. This transition is, if continuous, described 
by the Ising universality class. 

The observations of Chen et al>il are consistent with 
our analysis. At high temperatures, they observe a 
Coulomb phase (corresponding to the quantum super- 
fluid), as in the original model with full cubic symmetry. 
As the temperature is lowered, the system enters a 'para- 
magnetic' phase (the valence-bond Mott insulator), with 
no broken symmetry and confined monomers, via a tran- 
sition in the inverted- XY class. Symmetry between the 
x and y directions (Xi on the kagome lattice) is broken 
at a lower-temperature transition into the dimer crystal 
(quantum density- wave order), which is in fact found to 
be strongly first-order. By considering the interactions 
that can be added to the Lagrangian Co [Eq. (|47f ] in 
this reduced-symmetry case, Chen et alJ^. have shown 
that this first-order transition can be naturally under- 
stood as a "spin flop" , a reorientation of the vector order 
parameter <p. 
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